Exciton spin relaxation in single semiconductor quantum dots 
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We study the relaxation of the exciton spin (longitudinal relaxation time Ti) in single asymmetrical 
quantum dots due to an interplay of the short-range exchange interaction and acoustic phonon 
deformation. The calculated relaxation rates are found to depend strongly on the dot size, magnetic 
field and temperature. For typical quantum dots and temperatures below 100 K, the zero-magnetic 
field relaxation times are long compared to the exciton lifetime, yet they are strongly reduced 
in high magnetic fields. We discuss explicitly quantum dots based on (In,Ga)As and (Cd,Zn)Se 
semiconductor compounds. 

PACS numbers: 78.67.Hc, 72.25.Rb, 63.20.Ls 



The current interest in the manipulation of spins in 
semiconductors is based on the abiUty to control and 
maintain spin coherence over practical length and time 
scales. Of particular interest for possible applications 
in quantum computing is the storage of spins in zero- 
dimensional semiconductor structures. Recent theoreti- 
cal investigations of the single carrier-spin relaxation in 
semiconductor quantum dots (QDs) found that the elec- 
tron spin-flip transitions are very slow with relaxation 
times of milliseconds and longer P, Q. But since most 
concepts for quantum information processing involving 
QDs are based on optical generation, manipulation and 
read-out of spins, one has to investigate the dynamics 
of the excitonic spin. In our previous communication!^ 
(motivated by experimental results for InAs QDs0| at 
high temperatures T > 40 K), we studied the exciton- 
spin relaxation in QDs related to a second-order process 
governed by optical phonons. Further detailed theoreti- 
cal studies of the exciton-spin relaxation in QDs are still 
missing. 

Experimental investigations on the exciton-spin dy- 
namics in QDs refer mostly to the spin-coherence 
problem [3, la, i-^-, they determine the transverse relax- 
ation time T2. Studies of the longitudinal spin-relaxation 
time Ti are rare since they require strict resonant ex- 
citation conditions and/or high magnetic fields. Re- 
cent resonant-excitation experiments on InAs QDs indi- 
cate no exciton-spin relaxation at low temperatures: the 
exciton-spin is totally frozen during the radiative life- 
time even for high magnetic fields up to 8 tQ. Similar 
results are reported for CdSe QDs with large lateral di- 
mensions (comparable to twice the exciton Bohr radius of 
^ 5 nm in bulk CdSe crystals) |^ . However, longitudi- 
nal spin-relaxation times comparable to the exciton life- 
time were deduced from high magnetic field experiments 
for much smaller CdSe QDs[l3|. The observed magnetic 
field dependence (the relaxation rate strongly increases 
with magnetic fields) points on an acoustic-phonon me- 
diated exciton spin-flip as underlying mechanism. The 
same conclusion is drawn from the linear temperature 



dependence of the spin-relaxation rate reported for non- 
resonant excitation in InGaAs disks [13 ■ 

We propose in this paper an intrinsic mechanism for ex- 
citonic spin-flip transitions at low temperatures resulting 
from a deformation-induced exchange interaction. Our 
model calculations are able to qualitatively reproduce the 
above mentioned trends of the experiments. We study 
the dependence of the exciton spin relaxation time Ti on 
temperature, magnetic field, and the QD size. Finally, we 
draw conclusions on the suitability of various materials 
for the storage of optically generated spins. 

Extensive experimental studies have identified the 
main features of the exciton fine structure in self- 
organized QDs by means of single-dot spectroscopy 
Il3l |. Such QDs are usually strained and have an asym- 
metrical shape with a height smaller than the base size. 
It has been shown that a reduction of the QD symme- 
try lifts degeneracies among the exciton states and re- 
sults, in particular, in a splitting of the exciton ground 
state. Thus, as a consequence of strain and confinement 
in the growth direction of a QD, the ground states of 
the heavy-hole (hh) and light-hole (Ih) excitons are well- 
separated and the hh-exciton has the lowest energy. The 
hh- and Ih-exciton quartets are characterized by the pro- 
jections Jz = ±1, ±2 and — ±1, of the total an- 
gular momenta J = 1,2, respectively. The short-range 
exchange interaction splits the ground states of both hh- 
and Ih-excitons into doublet states (so-called singlet- 
triplet splitting), as is shown in Fig. ^a). The lateral 
anisotropy of a QD leads to a further splitting of the op- 
ticaly allowed doublets | ± 1) into two levels (labeled \X) 
and \Y)) with dipole moments along the two nonequiva- 
lent in-plane QD axes. This splitting is due mostly to the 
long-range exchange interaction and originates from the 
lateral elongation of the QD's0, 113. n^ . Continuous 
wave single-dot spectroscopy experiments have clearly 
evidenced the two related, linearly polarized optical tran- 
sitions e.g., in GaAs interfacial dots or self-organized 
InGaAs QD's0. Measured magnitudes of this splitting 
reach some tens or even hundreds of MeV|l|il[l|. Re- 
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FIG. 1: Schematic diagram showing (a) the sublevels of the 
heavy-hole-exciton and light-hole-exciton ground states in 
symmetrical QDs, (b) the phonon-assisted spin-flip processes 
within the radiative doublet of the heavy-hole-exciton in an 
asymmetrical QD, and (c) the Zeeman radiative doublet of 
the heavy-hole-exciton. Optically inactive states are shown 
by dashed lines. 



laxation processes between the radiative states changes 
the occupation of the exciton levels and are manifested 
in a change of the linear degree of polarization of the 
luminescence 0, 0] ■ 

In what follows we consider (longitudinal) spin- 
relaxation processes between the allowed \X)- and 
states of the hh-exciton. In such processes, the hh- 
exciton stays in the same spatial state and just flips its 
spin (i.e., the electron and the hole spin flip simulta- 
neously), by emitting or absorbing an acoustic phonon. 
Since the bare electron/hole -phonon interaction does not 
contain spin operators, it cannot directly couple the in- 
volved exciton states, however, this can occur via the hh 
and Ih exciton mixing due to the interplay of the the 
short-range exchange interaction and the lattice defor- 
mations. 

The short-range (isotropic) exchange interaction is 



given by the Hamiltonian 
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where (? and J are the electron spin and the hole to- 
tal angular momentum operator, respectively, and Ajf is 
the (singlet-triplet) exchange energy. Phonon-induced 
deformations come into play via the off-diagonal terms 
in the Bir-Pikus Hamiltonian IITI 
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where b and d are the exciton deformation potentials, 
Eij is the deformation tensor, and [JiJj] = JiJj + JjJi. 
(An overall shift of —ae of the Ih and hh exciton levels is 
omitted). 

In the basis of the hh- and Ih-exciton states {\ — l)hh, 
\ + l)hh and \ -l)ih, | + respectively), the total Hamil- 
tonian H = TCie) + Hex is represented by the following 
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The origin of energy axis is flxed at the hh-exciton states 
which are separated from the Ih- exciton by Eih, see 
Fig- da). The deformation-dependent off-diagonal term 
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anisotropic part of the splitting of the allowed doublets, 
since the energy difference Eih is, typically, much larger 
(of the order of several tens of a meV) than the flne struc- 
ture energies For the same reason, the hh-lh-exciton 
mixing can safely be regarded as a perturbation. It fol- 
lows from Eq. ^ that perturbations which mix hh- and 
Ih-exciton states result in a coupling of the | — l)hh~ and 
I + l)/t;i "States and, consequently, lead to the mixing of 
the 1^)- and \Y)- states of the hh-exciton. For QDs 
which are elongated in the [110] direction, this coupling 
is due to vibrations causing {sxx — £yy) deformations 
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Thus, a coupling of the \X)- and |y)-states is a result 
of an interplay between the short-range exchange inter- 
action and deformations: neither the exchange interac- 
tion nor the deformation perturbation j alone can couple 
these states. According to Eq. the matrix element 
for the relaxation processes among the \X)~ and \Y)- 
states of the hh-exciton is 



MxY = ^ b {{Exx - £yy)), 
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where the expectation value is taken with respect to the 
spatial part of the exciton wave function. The elastic 
strain is given by 
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where the phonon displacement field is 
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Here q , e and ujq — sq are the phonon momentum {q = 
\q\), polarization vector and frequency, respectively, V is 
the normalization volume, p is the mass density, and s is 
the sound velocity. Upon substitutions given by Eqs. ® 
and |(7J) we obtain for the matrix element in Eq. 
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where indices e and h stand for electrons and holes, Mor 
is the orbital part of the matrix element, and <i>o(r^,rX) 
is the exciton ground state envelope wave function. 

The spin-relaxation rates accompanied by phonon 
emission and absorption are given by Fermi's golden rule 
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where hft is the energy splitting between the \X)- and 
|y)-states (see Fig.HIb)), and N^^ ^ i/^^^hujjksT _ 
is the thermal phonon distribution function. 

We discuss now what determines the value of the or- 
bital part of the matrix element Eq. 0. Because of en- 
ergy conservation, \Mor\ = \M^^ + M^^] has to be cal- 
culated for the phonon momentum q — ujq/s = Q/s, 
while the characteristic extent of the orbital wave func- 
tion <i>o is, roughly speaking, the lateral size L of the 
QD. (For flat QDs, the QD hight < L.) Due to 
the oscillatory behavior of the exponential functions in 
Mor, we have IM^;!"] « 1, if f7L/s < 1. To illustrate 
this approximation, we consider the model of a strongly 
confined QD with the shape of an parallelepiped and in- 
finite barrier. Then the exciton states are reasonably 
well approximated by noninteracting electron-hole pair 
states. In addition, the single carrier-phonon matrix el- 
ement \M^^{q)\ = \M^^{q)\ = \Mor{q)\ separates in the 
x,y,z coordinates. For a confined direction (say, the x- 
direction), the following expression holds 
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where Qx — qxLx/2 and is the QD size in the x- 
direction. As seen from Eq. Hll|) . the matrix element 
|-^or(9a;)| IS Smaller than unity for any q^ ^ 0, but 
it approaches unity for small qxL^ <C 1 ( \Moj.iqx)\ ^ 
|1 — Ql/n'^l). With increasing qxLx > 1, \Mor{qx)\ de- 
creases rapidly ( \Mor{qx)\ ~ \ sinQx/Qx\)- Thus, if we 
replace the orbital matrix element in Eq. (jSJ by unity, 
\Mor{q)\'^ = 1, the exciton-spin relaxation rates Eq. ifTHI 
will be overestimated. Note that for typical values of 
and s we have (s/fi) ~ 40 nm, L is usually about 
10-20 nm^^, so that the above mentioned assumption 
Q,L/ s ^ 1 is reasonable. As a result, a lower limit for the 
relaxation times is given by 
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In order to obtain a numerical estimate, we use typi- 
cal parameters together with an estimated hh-lh- 
exciton sphtting oiEih = 10 meV, i.e., Eih/Agt — 50. As 
a result, we get for the spin relaxation time scale Ti = 
(t^"^ + ■^f^)"^ ~ ^^^^ T = 10 K) and Ti - 150 ns 



(at T = 100 K). Thus, for typical quantum dots and low 
temperatures, the zero-magnetic field relaxation times are 
very long compared to the exciton lifetimes of 1 ns. 
This is in qualitative agreement with resonant-excitation 
experiments at low temperatures 4, 8] . Please note, that 
because of the small inter level splitting Tifl ~ 0.1 meV, 
the relaxation times T| and are already comparable at 
T > 2 K. For higher temperatures Ti is inversely pro- 
portional to the temperature. But for T > 100 K, the 
interaction with LO phonons via a second-order process 
will be the dominant relaxation mechanism0, 0| . 

The calculated spin relaxation times according to 
Eqs. H12I13II are large partly because of the small phonon 
density of the states (which is ~ (ri/s)^) at the scale of 
the interlevel splitting hft. But, considerably large in- 
ter level splittings of hilB — QP'bB ~ 1.5 meV (where 
g is the exciton effective ^-factor) arise in high-field 
magneto-optical experiments. Here an enhancement of 
spin-relaxation rate is expected. The spin fiip in Fara- 
day geometry {B \\ z) occurs between the | — 1)4/1- and 
| + l)/ih" states of the hh-exciton (see Fig.mc))!;13|. Con- 
sequently, exciton-spin relaxation processes between the 
Zeeman sublevels are reflected in a change of the circular 
degree of polarization of the luminescence. 

For magnetic fields and zero temperature, a lower limit 
of the spin relaxation time is given by 
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The relaxation time in Eq. Hlt)|) varies as the square of 
the ratio Eih/Ast and depends strongly on the QD size 
because both the splitt ing E/ h and the exchange energy 
Ast are size-dependent pol l21| . If the splitting Eih is due 
mostly to the strain (large lattice mismatch) , Ti changes 
with the lateral and vertical size as and L^, respec- 
tively. For small lattice mismatch, Eih is due mostly to 
the vertical confinement and, therefore, Ti ~ L'^Lj^. In 
any case, the relaxation time Ti strongly decreases when 
decreasing the lateral size L of a QD. This result can ex- 
plain qualitatively the observations for the CdSe QDs of 
different lateral sizes in high magnetic fields. As noted 
above, no spin relaxation was found for large CdSe QDs 
(with L about 10 nm)Q, while for small CdSe QDs (with 
L about 3 nm - 5 nm) an efficient exciton-spin relaxation 
was observed 

To obtain numerical estimates for the spin relaxation 
in magnetic fields, we use the experimental values for 
Eih,Ast and g, and we take the typical values [l8l| for 
other factors in Eq. (I14|. For InAs/GaAs QDs, we use 
Ast = 0.2 meV, g = 30 and Eih = 30 meV {Eih of 
several tens of a meV is reported in Ref . [l^ ) . As a result 
we get a lower limit for the relaxation time Ti ~ 13/i?'^ 
(/iS T'^) which is much larger than the exciton lifetime 
scale '--^ 1 ns even at high B (e.g., Ti ~ 25 ns at B = 8 T). 
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Indeed, no spin relaxation is found for fow temperature 
experiments in high magnetic fields (Jj. 

This case is different for small CdSe/ZnSe QDs investi- 
gated in Ref. . Here, values of Ast ~ 1 meV, g ~ 2 
and Eiii ^ 40 meV are reported. As a result we calculate 
Ti ~ 5 ns at i? = 6 T which is close to values of 2.5 
ns - 3.4 ns deduced from the experiments [Td|. For large 
CdSe/ZnSe QDs investigated in Ref. 9J no data on the ex- 
change energy Ast are reported but we still can estimate 
the relaxation rate. The lateral size of the QDs in this 
case is about a factor of 2.5 larger than in The con- 
sequence is a much smaller exchange splitting Ast ~ 0.2 
meV |2ll | which results in much longer relaxation times 
of ~ 100 ns (at i? = 6 T). As already mentioned, no spin 
relaxation is observed in the experiments of Ref.0. 

Let us now consider the magnetic field dependence of 
the spin-relaxation time which was measured in small 
CdSe/ZnSe QDs [13 • The experimental results indicate 
a roughly anti-proportional decrease Ti oc up to 

8 T. As seen from Eqs. (|14I16|I . Ti scales as B^^ for 
small B (r^B < Q), in agreement with the experiment. 
But it would depend much stronger on the magnetic field 
(as B^^) for large B. This difference to the experiment 
can originate from the treatment of the orbital matrix 
element, which up to now has been fixed as 1. But a 
factor of l/|Mor-p has to be accounted for explicitly in 
Eq. (|14|l for high magnetic fields. The reason is, that the 
relevant phonon momentum q{B) ~ fis/s is linearly pro- 
portional to B. Thus, l/|MorP increases with magnetic 
field and becomes important at q{B)L > 1. This is the 
case e.g., for B > 5 T for CdSe QDs with a lateral size of 
L ~ 4 nm and is the regime of the experiments in lO] . So 
we can confirm the general trend of the experiment _^10J, 
namely that Ti decreases with magnetic field. An analyt- 
ical dependence of Ti on B (at high B) would require an 
explicite calculation of |MorP with the QD exact param- 
eters, which for the QDs used in Ref. 10], unfortunately, 
are not known. 

There are many experiments, which determine the po- 
larization degree of the luminescence from QDs under 
nonresonant excitation conditions. We just want to dis- 
cuss the results of Ref.^| as a typical representative. 
Rather short polarization decay times on the order of 1 
ns have been found here. But definite conclusions on the 
underlying spin-relaxation mechanisms are not possible 
in a scenario where the electron-hole pairs relax through 
a multitude of barrier states and/or excited quantum-dot 
states. To identify spin-relaxation processes, the split- 
ting of the exciton ground state has to be large (as for 
large magnetic fields) or strictly resonant excitation has 
to be used (see U]). 

The above discussions show that a direct quantitative 
comparison of our proposed relaxation model to exper- 
imental data is at the moment difficult at best. One 
obvious reason is that the relaxation times, once they 
are much larger than the exciton lifetime, cannot be 
quantified experimentally. The second reason is, that 
only very few experiments under strict resonant excita- 



tion and/or in high magnetic fields have been performed. 
But still, there are a couple of experimental trends which 
are explained consistently by our model. The value and 
magnetic-field dependence of Ti in small CdSe QDs can 
be reproduced. We correctly predict the experimental 
findings that long relaxation times are expected for InAs 
QDs as well as large CdSe QDs at low temperatures even 
in strong magnetic fields. 

The experimental differences found for InAs vs CdSe 
QDs and for large vs small CdSe QDs strongly support 
our proposed model. Let's consider alternative spin- 
relaxation mechanisms within the radiative doublet. Pro- 
cesses relying on the exciton motion like for quantum- 
well excitons 0] are suppressed in QDs. Assume, the 
electron and hole would individually flip their spins. In 
contrast to our proposed mechanism, this is a higher- 
order process involving intermediate dark exciton states 
(the excitons finally have to return to radiative states 
to contribute to a depolarized luminescence). The re- 
laxation time here is determined by the longest spin- 
flip time for a single carrier (electron or hole). The spin 
flip of the electron has to rely on mechanisms based on 
the spin-orbit coupling. These relativistic effects have 
been shown by Khaetskii and Nazarov [ll| to be inef- 
ficient in QDs. Finally, the spin-orbit interaction and 
carrier-phonon coupling are not different enough in InAs 
and CdSe or for large and small QDs to explain that 
the spin relaxation times change by orders of magnitude. 
This dramatic change found in the experiments is, how- 
ever, readily explained by the strongly differing exchange 
splitting in InAs and CdSe and its significant enhance- 
ment in small dots. Since the exchange splitting enters 
quadratically in our proposed relaxation scheme, we are 
quite convinced to have identified the most prominent 
spin-relaxation process for QDs. 

Finally, we want to discuss different material systems 
with respect to their suitability for spin storage. The rel- 
evant parameters for several semiconductor compounds 
forming QD's and their matrices are presented in TableQ] 
Among the different material parameters in Eq. (|15|l . the 
characteristic energies Eih and Ast are strongly affected 
by the structural composition. Thus we can only give 
general trends. We will consider QDs with typical sizes 
of about Lz ~ 2-5 nm and L ^ 10-20 nm which are 
in the strong vertical confinement regime. Moreover, 
InxGai_xAs QDs are also strongly confined in the lat- 
eral plane, and the exchange splitting Agt in these QDs 
significantly exceeds its bulk value of A^j 21]. In CdSe 
QDs the bulk value of the exchange splitting is already 
quite large. In CdSe QDs with lateral sizes larger than 
the exciton radius of ~ 5 nm the exchange splitting is 
additionally affected by the vertical confinement. Since 
the exchange interaction in (In,Ga)As crystals is much 
weaker than in CdSe crystals, Agt is of the same order of 
magnitude p3] in strongly confined luxGai-xAs QDs and 
large, weakly confined CdSe QDs. Considering addition- 
ally the strain-induced Ih-hh splitting Eih, which has to 
be large to suppress spin relaxation, InAs QDs have the 
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TABLE I: Semiconductor compounds forming QDs (first column) and matrices (sixth column), A^t and aex axe the exchange 
splitting and the exciton Bohr radius in bulk compounds forming QDs, A^j gives typical values of the exchange splitting in 
QDs with sizes of Lz = 3 nm and L = 15 nm. mo/m = mo{l/mih — 1/mhh), where mih and ruhh are the Ih- and hh-masses, 
respectively. |Aa|/a = (a^ — am) /0.5{ad + am) and AEg = Egm — Egd with lattice constants and am and energy gaps Egd 
and Egm in the QD and matrix materials, respectively. Material constants are taken from Ref . except A^^ and a^^ for InAs 
crystals which are taken from Ref. 25. For the case of CdSe crystals, A^t, a^x, AEg, and mo/m are given for hexagonal CdSe 
crystals. 
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AEg 
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(nm) 
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(eV) 


InAs 


38 


0.3 


200 


37.5 


GaAs 


0.069 


1.1 


InAs 


38 


0.3 


200 


37.5 


Alo.3Gao.7As 


0.068 


1.8 


GaAs 


11 


20 


250 


9.3 


Alo.3Gao.7As 


0.0005 


0.48 


Ino.6Gao.4As 


24 


8.2 


220 


17.7 


GaAs 


0.042 


0.62 


Ino.4Gao.6As 


22 


12.1 


230 


13.8 


Alo.5Gao.5As 


0.027 


1.22 


Ino.7Gaa.3As 


31 


6.2 


210 


20.6 


GaAs 


0.049 


0.73 


CdSc 


5.4 


130 


450 


4.7 


ZnSe 


0.066 


0.9 



longest relaxation time Ti. 

Less suitable for spin storage is the GaAs/Alo.sGao.yAs 
system since its lattice mismatch is very small. Thus, 
among the QD structures listed in Table I, the shortest 
Ti relaxation times can be expected for strongly confined 
GaAs/Alo.sGao.rAs QDs. Likewise, small Ti times are 
predicted for small CdSe QDs since in the strong lateral 
confinement regime, A^t is considerably enhanced. 

In conclusion, the exciton-spin relaxation within the 
radiative doublet of the exciton ground state in single 
asymmetrical QDs is studied. As a possible intrinsic 
mechanism for such a process, the exciton spin-acoustic 
phonon coupling via the strain-dependent short range ex- 
change interaction is proposed. For zero-magnetic fields 
and low temperatures, the calculated Ti-relaxation times 
for typical QDs are long compared to the exciton life- 
time. A strong reduction of the relaxation times occurs 



in QDs in high magnetic fields. Nevertheless, numerical 
estimates for InAs QD's and large CdSe QDs give large 
values for the Ti -relaxation times even for high B ^ 8 T 
(up to a few of tens of nanoseconds). In addition, the re- 
laxation time Ti strongly decreases in strongly confined 
(in the lateral plane) QDs. The Ti -relaxation times esti- 
mated for small CdSe QDs in high magnetic fields reduce 
to the nanosecond scales. These predictions of our model 
are in qualitative agreement with experimental findings. 
Within the considered mechanism, we conclude that InAs 
QDs and CdSc QDs with typical lateral sizes of i ^ 10 
- 20 nm display a weak exciton spin relaxation, whereas 
for typical GaAs/Alo.sGao.rAs QDs and small CdSe QDs 
the Ti times are expected to be rather short. 

This work was supported by the Center for Functional 
Nanostructures (CFN) of the Deutsche Forschungsge- 
meinschaft (DFG) within project A2. 
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